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We investigate the dynamical properties of nanolasers comprising a few two-level emitters coupled 
to an optical cavity. A set of rate equations is derived, which agree very well with a solution of the 
full master equation model and makes it simple to investigate the properties of the system. Using a 
linearized version of these rate equations, we can analytically express the response of the nanolaser 
to a modulation of the pumping rate. These results are compared to the modulation response 
obtained directly from the master equation using a novel method. Using the rate equation method, 
we calculate the modulation bandwidth and show that, contrary to conventional semiconductor 
lasers, the nanolaser is typically over-damped and displays a dip in the modulation bandwidth as 
the two-level systems become inverted. Both these features can be traced back to the modeling of 
the emitters as two-level systems that are incoherently pumped. 
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The laser has evolved from table size apparatuses to 
truly nano sized devices, in much the same way that 
computer chips have been continuously minimized. The 
single or few-emitter nanolaser represents an extreme in 
terms of size. The emitters can be either atoms or quan- 
tum dots (QDs) that are coupled to a high-Q optical cav- 
ity, with mode volume on the order of the wavelength of 
light cubed. In these systems, the low number of emit- 
ters is compensated by a strong light-matter coupling 
and high Q values. Indeed, these nanolasers operate in a 
regime where experiments in cavity quantum electrody- 
namics (cQED) previously have shown strong-coupling 
behavior such as Rabi-splitting [TH3]. 

The experimental realization of such nanolasers is rela- 
tively new. For atoms the experiment may be performed 
by letting the atoms fall through the cavity and thereby 
obtaining a strong light-matter coupling during the tran- 
sit period @]. For lasers based on semiconductor QDs, 
the optical cavity is in the form of a planar photonic 
crystal [6] or a micropillar [7] . A strong coupling be- 
tween the QDs and cavity is obtained either by defining 
the cavity around a preselected QD [51 , or by select- 
ing a device with good spatial and spectral match from 
many samples. Although such lasers are in their infancy, 
they may have a promising future in applications where 
their small size and low power consumption are impor- 
tant, e.g. for on-chip optical interconnects and photonic 
integration. Other examples of lasers with one or a few 
emitters include ions in optical traps [5] and artificial 
atoms in quantum circuits |10j . 

In the simplest model for such a nanolaser, the emit- 
ters are modeled as two-level systems, which are coupled 
to a single mode of the electromagnetic field in an optical 
cavity. The two-level emitters are pumped incoherently, 
i.e. by using a reverse decay rate from the ground to the 
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excited state. Such models have previously been investi- 
gated [TTH2"U] and enjoy popularity, especially because of 
their simplicity. 

The optical cavity changes the local density of op- 
tical states experienced by the emitters. This means 
that light-matter interactions are suppressed for some 
frequencies and Purcell enhanced [3TJ [55] for others, 
namely for those at the cavity resonance. Such lasers may 
become thresholdless, meaning that upon pumping the 
laser will start lasing, without the abrupt change in out- 
put power and coherence known from conventional lasers 
[20l [23] . However, for an incoherently pumped two-level 
emitter, a large pumping rate will dephase the transi- 
tion to such an extent that the output will eventually be 
quenched and the coherence will be lost [TTJ [T5J. Thus, 
the output of a few-emitter nanolaser can be divided into 
three regimes: the low pump regime where the emission 
is either anti-bunched or strongly bunched depending on 
the number of emitter, the laser regime at intermediate 
pumping rates where the emission has a poissonian dis- 
tribution, and a quenched regime with chaotic emission 
at larger pumping rates. 

In this paper we derive a rate equation model for 
a nanolaser system with an arbitrary number of single 
emitters coupled to a cavity. This derivation is based on 
the Von Neumann master equation with Lindblad terms 
|24) . which is widely used in cQED models. The rate 
equations are easy to solve and give physical insight, as 
opposed to numerical solutions of the full master equa- 
tion. We proceed by using these rate equations to analyt- 
ically investigate the response of a nanolaser when sub- 
ject to a modulation of the pumping rate, and thereby 
gain a deeper insight into the fundamental dynamics of 
few-emitter nanolasers. Modulation responses have pre- 
viously been investigated both experimentally [25] 
and theoretically 27, 28J for lasers containing many emit- 
ters, but the present paper is the first to calculate these 
properties for system with only a few emitters. Previous 
investigations have predicted that nanolasers can exhibit 
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FIG. 1: Schematic of the system considered where N emitters are 
coupled to a single optical cavity with the coupling strengths gi. 
The cavity decays at the rate k and the emitters are subject to a 
background decay rate 7; and an incoherent pumping rate Pi. 



a very large modulation bandwidth caused by the Purcell 
enhanced spontaneous emission rate [SJJ HE] • This effect 
should be even more pronounced in few emitter lasers, 
since they rely on a large Purcell enhancement to obtain 
lasing. We calculate the modulation bandwidth for the 
nanolaser and show that it displays a characteristic dip 
when the system is lasing. We also show that the re- 
sponse of the system to a modulation is over-damped, 
whereas conventional semiconductor lasers typically are 
under-damped. The models developed in this paper can 
be used to assess the applicability of nanolasers in on-chip 
communication systems. 



I. MODEL 

A two-level model with incoherent pumping is one of 
the simplest possible models for quantum emitters and 
has enjoyed large popularity for analysing cQED systems. 
It has therefore been used extensively to model single 
emitter lasers [TTHT31 fTMTTl [151 |2"5] . Incoherent pump- 
ing models the excitation of an electron from the ground 
to the excited state, as e.g. from one atomic shell to an- 
other. In a semiconductor perspective this can be seen 
as a correlated creation of a hole in the valence band and 
an electron in the conduction band, i.e. the creation of 
an exciton. Although this model has also been used pre- 
viously for modeling laser based on semiconductor QDs 
[TBI HI"] , it lacks several features which are relevant in 
QDs. For a more realistic model of QDs, it may be im- 
portant to include additional levels with a more sophis- 
ticated pumping method and effects such as Coulomb 
interaction [30 33 and hybridization with states of the 
wetting layer [Ml [35]. However, a two-level model can 
be used to model atomic system which are generally free 
of the environmental effect of semiconductors, and may 
also be used to model lasing in systems consisting of su- 
perconducting qubits [ini [35] ■ Nevertheless, the general 



features that we predict and analyze should be relevant 
also for QD based lasers. 

The nanolasers considered here contain N emitters, 
modeled as two-level systems and coupled to a single 
mode of a cavity. See Fig. [TJfor a schematic. The light- 
matter coupling strength between the i'th emitter and 
the single mode of the optical cavity is gi, and the emit- 
ters are pumped individually with an incoherent pump- 
ing term with the rate Pi. Additionally, the emitters 
can decay with a rate 7,, which models spontaneous de- 
cay to other modes than the primary mode of the cavity, 
and also accounts for any non-radiative relaxation of the 
emitter. Finally, light is coupled out of the cavity at the 
rate k. 

To model the different decay channels, the system 
is coupled to a number of reservoirs. These reservoirs 
are assumed to be time-independent and memory-less. 
In this case the Born-Markov approximation can be 
made and the resulting master equation for the combined 
emitter-cavity system is [24] 



dp(t) 
dt 
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where p(t) is the density matrix. The Hamiltonian is of 
the Jaynes-Cummings form 

N N 

Hj C = twea' fru i al<Ti + y^ j hgi (a^ + ajaj , (2) 

i=i i=i 

where a\ (er^) is the raising (lowering) operator for the 
i'th emitter, which has frequency uj i: and where a and a* 
are the creation and annihilation operators for the sin- 
gle photonic mode of the cavity with frequency u> c . We 
shall assume uji = uj c , corresponding to zero detuning be- 
tween the emitters and the cavity. For identical emitters 
it is possible to tune the emitters and cavity into reso- 
nance. In the case of inhomogeneous broadening, such as 
for Stranski-Krastanov grown QDs, the cavity may act 
to select those QDs that are in resonance, but a finite 
detuning may need to be included. 

The Lindblad terms in the master equation are 

C-Kp{t) = I (a)ap{t) + p{t)a)a - 2ap{t)(J) (3) 



for the cavity decay and 



C x p{t) = 



N 

£ 
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-2F«p(t)r 



(4) 



where Cpp(t), with Tp = crt, represents the pumping 
process where the emitter is excited from the ground to 
the excited state. The background decay, which includes 
spontaneous emission to other modes than the cavity 
mode and non-radiative recombination, is represented by 
£ 7 p(i) with ri^ = Cj. 
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FIG. 2: (color online) a) The mean cavity population as a func- 
tion of the pumping rate for systems with between one and four 
emitters. The full lines are obtained using the rate equation model 
and the marks are numerical solutions to the master equation ob- 
tained using direct time-integration, except for the highest (marked 
with filled-in squares), which have been obtained using a Quantum 
Monte Carlo method, b) The calculated second order correlation 
function, </ 2 )(0) for the four cases. Three regimes are marked on 
the abscissa axis of b), corresponding to I: strongly bunched/anti- 
bunched emission, II: lasing and III: chaotic emission. The pa- 
rameters are g% = g = 300 ns" 1 for the light-matter coupling, 
7i = 7 = 10 ns — 1 for the background decay rate and k = 100 ns _1 
for the cavity decay rate. 



From Eq. (JlJ it is easily shown that the time evolution 
of the expectation value of a general operator A is [TB] 



d(A) d . A , 



= l<[Affjc]} + ?([aU] a + at [A, a]) 



r(<) 



(5) 



It shows that the time-evolution of an expectation value 
(A) is in general linked to the time-evolution of higher- 
order expectation values, such as e.g. (AcJa) and 
(atAa). This expression will be used to derive the rate 
equations in the next section. 

By using a basis consisting of the levels of the emitters 
and the number states of the cavity, the master equation 
can be transformed into a set of coupled linear differential 
equations. The Fock states of the cavity are truncated, so 
that only the lowest N c states are included. The choice 
of N c depends on, but is generally higher than, the mean 



population of the cavity. The resulting set of differen- 
tial equations may then be solved using direct integra- 
tion. However, when including several emitters and many 
states of the cavity, the Hilbert space grows rapidly, mak- 
ing direct integration unfeasible due to memory require- 
ments. In this situation a Monte Carlo method [37l [38] 
may be employed, which trades run time for greatly re- 
duced memory usage. For the Monte Carlo calculation 
we have used the software package QuTiP [3"§] . 

In Fig. [2] a) the mean cavity population of systems 
with between one and four emitters is shown for different 
values of the pumping rate P. The Monte Carlo method 
has been used for the system with four emitters and high 
cavity population, and these results are shown with filled- 
in squares. For low pumping rates the cavity is scarcely 
populated (marked as regime I in the figure), but as the 
pumping rate is increased the mean number of photons in 
the cavity grows proportionally to the pumping rate P. 
This increase in the population of the cavity brings the 
system into the lasing regime, where the cavity photons 
stimulate the emission of further photons 11 (regime II 
in the figure). Since the systems are nearly thresholdless, 
it is hard to distinguish regime I and II using the cavity 
populations alone. However, as will be discussed later, 
the statistics of the cavity field can be used to identify 
the regimes. For very large pumping rates the output 
quenches, i.e. the number of photons decreases towards 
zero, because the incoherent pumping dephases the tran- 
sition from the emitter to the cavity (regime III in the 
figure) . This phenomenon has been described before and 
is referred to as self-quenching [TTJ [T2J Q13 [T7]. These 
three regimes are general for systems with any number 
of emitters, as illustrated in Fig. [2] where the cavity pop- 
ulations for up to four emitters are shown. When more 
emitters are coupled to the cavity, the cavity population 
can grow larger before the quenching sets in. 

Reaching the lasing regime is conditioned on having a 
good cavity with a low cavity decay rate, ft, and a large 
light-matter coupling, g. For a system with N identi- 
cal emitters the collective light-matter coupling strength, 
9 n = gVN, increases with the number of emitters |40) . 
Therefore, a system which cannot reach the lasing regime 
with a single emitter may be able to do so if more emit- 
ters are coupled to the cavity. 

Laser light is conventionally associated with a coher- 
ent state of the optical field, however it is well known 
that a coherent state is difficult, if not impossible, to ob- 
tain in quantum optics |41| . It can be shown that the 
cavity state of the systems described here can never be 
fully coherent by writing the equation of motion for (a) 
using Eq. ([5]). This shows that (a) depends only on 
(a) which in turn couples to higher order operator av- 
erages such as (aa'aj. But there are no direct driving 
terms for (a), (a) or any of the other operator averages 
in the hierarchy, which means that in steady state (a) 
is zero. A coher ent state, on the other hand, always has 
I (a) I = \J (at a). This conclusion has been reached previ- 
ously by examining the structure of the master equation 



4 



in Eq. ([IJ, where it can be seen that there is no cou- 
pling between different Fock state of the cavity when the 
emitters states are identical [2]. This means that the 
reduced density matrix for the cavity states, obtained by 
tracing out the degrees of freedom associated with the 
emitters, p cav = Tr cm i t [p], only has diagonal elements. 
This is different from the density matrix of a coherent 
state H21 



— \/n\m\ 

n , m v 



(6) 



which, has non-zero off-diagonal elements. However, the 
system can be in a state where the diagonal terms of the 
density matrix are similar to those of a coherent state, 
but the others are zero [H] . Such a state is characterized 
by poissonian photon statistics. 

Formally the photon correlation functions of all orders 
are required to fully identify the photon statistics of the 
cavity field as poissonian, but the second order correla- 
tion is a good indicator of the statistics 43 . The sec- 
ond order correlation function is defined as, g^(0) — 
{(n 2 ) - (n)) / (n) 2 . When g (2) (0) < 1 it indicates that 

the emission is anti-bunched and when </ 2 )(0) = 2 it 
signifies chaotic light, similar to that emitted by an in- 
candescent light bulb. A poissonian photon distribution 
is characterized by [42] g( 2 '(0) = 1, which is thus the 
criterion we use for lasing. 

In Fig. [2] b) the second order correlation is shown, 
and its value has been used to identify the three dif- 
ferent regimes shown at the bottom of the figure. We 
see that in the low pump limit (regime I in the figure) 
the single emitter laser has a anti-bunched cavity field, 
and thus works as a single photon source. The cause of 
the anti-bunching of the emission is that the single emit- 
ter must be reexcited between two consecutive photon 
emissions, thus the probability of having a two-photon 
state in the cavity is small. Conversely, a system with 
more than one emitter can emit two photons into the 
cavity at the same time, and therefore shows a strong 
bunching behavior. We see that the system with three 
emitters shows a smaller bunching than the systems with 
two and four emitters. This dependency on the parity of 
the number of emitters is a sign of cooperativity, which 
means that the bunching is not accidental, but an indi- 
cation that the emitters interact through the cavity field 
[15] . As the pumping rate is increased the emission of 
all the systems tends towards </ 2 )(0) = 1 representing 
lasing (regime II in the figure), but without showing a 
clear kink in the input-output curves, which is normally 
used as the criterion for reaching the threshold. As de- 
scribed previously, this is well known for lasers where a 
large fraction of the spontaneous emission goes into the 
lasing mode [20] [23] . As the pumping rate is increased 
further the systems experience self-quenching (regime III 
in the figure), as described above, and the cavity field 
becomes chaotic with </ 2 )(0) = 2. 



II. THE RATE EQUATION APPROXIMATION 

Rate equation based models are very well estab- 
lished for modeling of conventional semiconductor lasers 
[HISS]. However, conventional rate equations cannot im- 
mediately be applied to few emitter nanolasers. Indeed, 
a single emitter that is strongly coupled to a high quality 
cavity, i.e. with similar parameters as required for las- 
ing, is known to operate in the regime of coherent Rabi 
oscillations. Such coherent oscillation usually cannot be 
predicted by a rate equation formalism. However, as we 
shall show, the large amount of dephasing which is in- 
duced by the incoherent pumping destroys the coherence 
and therefore enables a rate equation model description. 

In deriving a set of rate equations, we start by writing 
the equation for the mean cavity population using Eq. 
HJ) with A = a)a: 

|<ata> = |Tr[ a t ap ( t) ] 

N 

= 2 ft Im [<^ at )] - K < ata ) • ( 7 ) 



Similarly we write the excited state probability of the 
Vth emitter as 



<j\<Jip{t) 



--p, 



o-io-j) - 2gilm [(^a 1 )] - 7» (ofo) • (8) 



We see that both (a^a) and (^ajaij depend on the 

photon assisted polarization {<Ji(v\, which likewise can 
be derived from Eq. |5) as 



o]ajaa)) — (a; a] a* a 



(9) 



At this point we make two approximations. 1) 
we set 4z (o"ia^) = 0, which means that the pho- 
ton assisted polarization adiabatically follows the higher 

order operator averages (^x\aiaa)j and ^a i aja i a^j . 

This is a good approximation when the system is 
not operating in the coherent regime and the dephas- 
ing is large. 2) Using cluster expansion [3(31 [46], we 

set (^a\(Jiaa)^ = ^ejei^ (sa^) + (aiCi^) ^°l a ^ ana - 

(^Jia\a)aS = ^er,;crj^ (a^a) + (aicfl) (^y\a^. The approx- 
imation consists of ignoring the quadruplet correlations 
[3"U] 5 (ajaiaa^ and 5 (ai(j\a) a*J , which is justified when 

the dephasing of the photon assisted polarization, given 
by I (Pi + ji + k) , is large. This specific application of 
the cluster expansion method is known as the doublet 
approximation [3U1 |4"B1 ITT] . 



5 



By solving for (<7;a^) in Eq. (j9j) and inserting into Eq. 
([7]) and pL we arrive at the rate equations 



dn 
~dt 



N 

~ = E (^("^ - + - «n„ (10) 



dn, 



dt 



= P in ® - F,(n« - n«)n a - (F + 7i )n« (11) 



where rig = (<Ji<j\ 



and rie = ('cjcr, ) are the proba- 



bilities that the «'th emitter is in the ground or excited 
state respectively, and where n a = (ala) is the mean 
cavity population. Since we assume that the electron is 
confined to the two-level emitter, it must be either in the 
ground or the excited state. Thus n e and n g are related 

by n g i] + rii ] = 1. 

The factors Fj, which enter the rate equation are given 

by 



4.9? 



Pi 



■7i 



(12) 



and are the Purcell enhanced rates into the cavity. They 
specify the coupling between the cavity and the i'th emit- 
ter. Importantly, we notice that they depend on the 
pumping rate, which enter the expression in a similar 
way as the decay rates. 

The interpretation of the different terms in the rate 
equations is straight forward. The first term of Eq. (10 1, 



(i) 



)n ai corresponds to stimulated emission or 
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Fi(n 

absorption. If the emitter is not inverted, i.e. rie' < n g 
the term is negative which corresponds to absorption. 
Once the emitter is inverted the term becomes positive 

and corresponds to stimulated emission, which is respon- 

(i) 

sible for the laser action. The second term Fin e cor- 
responds to spontaneous emission, which like the stim- 
ulated rate has been enhanced by the presence of the 



cavity [22]. The term Pn g is the pumping term and 
depends on the occupation of the ground state. 

The self-quenching seen in Fig. [2]can now be explained 
by considering Fj. The denominator in Eq. (12 1 de- 



scribes the dephasing of the transition, which becomes 
large, and is responsible for the decrease of the cav- 
ity population, as the pumping rate is increased. The 
mean number of photons in the cavity is a result of two 
competing processes. On one hand a higher pumping 
rate increases the excited state population of the emit- 
ters thereby increasing the stimulated emission into the 
cavity. On the other hand the increased pumping rate 
give rises to dephasing which decreases the coupling to 
the cavity. As the emitters become fully inverted, the 
stimulated emission rate cannot be increased further and 
the output quenches. In the limit of very large pumping 
rates, we see from Eq. ( 12 ) that F — ► and the cavity 



population approaches zero [TTHT31 ITS] . 

The solution to the rate equations for one to five emit- 
ters is shown in Fig. [2] a). The results are compared with 



the values from the full model described in the previous 
section and, as can be seen, the two models agree well. 
In particular, the agreement is very good for large values 
of the pumping rate where, as discussed above, the large 
dephasing rate eliminates the higher order correlations 
between the emitter and the cavity, thus improving the 
approximations. 

We note that from Eq. Q and (|8| the following steady 
state relation can be derived 



^(Fi+7iK 



N 

( =£f. 



(13) 



This relation has previously been derived for the one 
emitter case (N = 1) [T31 US], and is a very powerful 
relation which we will use later. It should also be noted 
that this expression is derived directly from the master 
equation, without the approximations used in deriving 
the rate equations. 



III. MODULATION RESPONSE 

The modulation response of a nanolaser can be cal- 
culated in several ways. One can directly induce a small 
periodic modulation of the pumping rate when the laser is 
in steady state, and simulate the response in dependence 
of the frequency. Alternatively, the linear response can 
be calculated directly from the equations governing the 
system. In this section we will use the rate equation to 
calculate the linear response of the system, and thereby 
obtain an analytical expression for the modulation band- 
width. In appendix [S] we outline how to calculate the 
modulation response directly from the master equation 
and we use this method to validate the rate equation 
based method. 

The behavior of the system is generally non-linear, 
but its small-signal response can be analyzed from a lin- 
earized version of the rate equations. They are linearized 
around the steady state values n a and n e , by setting 
n a =n a + Sn a and rie = rie + 5n£ • The pumping term 
is expanded around the steady state value by writing it 
as Fj = Pi + SP, and the rate F i; which also depends on 



Fi + j^SP 



F, 



,5P. 



the pumping rate, as F — ^ , dp „^ — ^ , 4g2 
In matrix form the linearized rate equations for a single 
emitter (N — 1) become 



d 
dt 



5n a 

dn,. 



7ea lee 



5n a 
5n P 



5P 



-jp 



1p 



(14) 



This is a linear system driven by the small perturbation, 
5P, in the pumping rate. The element 7 ae = F(2n a + 1) 
determines the response of the cavity population to a 
variation of the emitter occupation. The first term orig- 
inates from stimulated emission/absorption and the sec- 
ond term is the contribution from spontaneous emission. 
The other off-diagonal element, 7 eQ = F(n e — fig), de- 
termines the response of the emitter caused by a small 
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perturbation in the cavity population, and originates 
from the stimulated emission/absorption term in the rate 
equations. The diagonal elements are i aa = j ea — k 
and 7 ee = — i ae — 7 — P and contain the same elements 
as the off-diagonal terms, in addition to terms originat- 
ing from the decay of the cavity and emitters. Finally 



contains the contribution to 



IP = f^r [(n e - n g )n a 
the variation of n a and n e originating from the fact that 
F varies with P. 

With the linearized version of the rate equations it is 
easy to derive the modulation response. By taking the 



Fourier transform on both sides of Eq. ( 14 ) and solving 



for Sn a (uj), the modulation transfer function, H(u>) 
^P^ttX, can be written as 

on a (0) 5 



H(u) 



^aa 1 at 
~^fea Tee 



-7p 



7<xe 



n g + 7P Tee + iu 



7 QQ + iuj j ae 

~lea lee + «W 



-IP lae 
rig + IP lee 



(15) 



In Fig. [3] the absolute square of the modulation trans- 
fer function \H(lu)\ 2 , referred to as the modulation re- 
sponse, is shown for three different values of the cavity 
decay rate. The modulation response obtained using Eq. 
( 15 1 is compared with numerical results obtained directly 
from the master equation, as outlined in appendix|Aj The 
rate equation model can be seen to agree quantitatively 
with the master equation results. 

From Fig. [3] it can also be seen that the modulation re- 
sponse decreases monotonically without the appearance 
of the usual relaxation oscillation resonance observed for 
conventional semiconductor lasers [H]. This can be in- 
vestigated by calculating the discriminant for the char- 
acteristic equation of the matrix in Eq. ( 14 1 which is 



D = {lae ~ leaf + 2 (T - «) (lae + lea) + (T - nf (16) 

where r = P + 7. If the discriminant is positive the sys- 
tem is over-damped, and the solutions are exponentially 
decaying functions, whereas if D is negative the systems 
is under-damped which means that the solutions are os- 
cillatory, and a resonance may be seen in the modulation 
response. Since j ae + j ea is positive, D can only be neg- 
ative if P < k — 7. Therefore, with the high pumping 
rate required for lasing, the system will always be over- 
damped while in the lasing regime. 

The discriminant for a conventional semiconductor 
laser is similar to Eq. (16 1, with the exception that T 
is independent of the pumping rate. It is instead given 
by 7s P which, similar to 7, is the spontaneous decay rate 
to other modes than the lasing mode and is independent 
of the pumping rate when the system is above the las- 
ing threshold [44] . Since 7 sp typically is much smaller 
than k, the discriminant can be negative, signifying an 
under-damped system. 

The main difference between the single emitter laser 
and the conventional semiconductor laser is thus the in- 
clusion of the incoherent pumping rate P in the second 
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FIG. 3: The modulation response |_ff(o;)| 2 for a single emitter as a 
function of the frequency f = lu/2tt for three different values of the 
cavity decay rate re. The modulation resp onse is calculated from 
the linearized rate equations in Eq. ( |15| (full line) and directly 
from the master equation using the method outlined in appendix 
[A] (marks). The horizontal line indicates a 3 dB decrease from the 
zero frequency response. The parameters are P = g = 300 ns _1 
and 7 = 10 ns _1 . 



term of the discriminant. This term stems from the way 
the pumping is represented in the model. In Eq. ( 11 ) the 



pumping rate is Pn g which means that the pumping is 
dependent on the state of the emitter. This is contrary 
to a conventional semiconductor laser where the pump- 
ing term is just P. Using the relation n g + n e = 1, we 
can write the pumping term as P(l — n e ), which gives the 
conventional population independent pumping rate, but 
also shows that there is an additional decay rate of the ex- 
cited state, —Pn e , as described above. This term enters 
in the same way as in e and therefore also gives rise to 
dephasing. The term is responsible for the P-dependcncc 
of r and is thus the origin of the fundamental difference 
in the dynamics between a conventional semiconductor 
laser and a nanolaser modeled as an incoherently pumped 
two-level system. 



IV. THE MODULATION BANDWIDTH 

Figure [3] shows that the intensity of the single emitter 
laser can follow the modulation until a certain frequency 
after which the response decreases drastically. The 3 dB 
frequency is the frequency at which the modulation re- 
sponse has decreased to half of its initial value and is 
also denoted the modulation bandwidth. In Fig. [3] we 
see that as k is increased the modulation bandwidth also 
increases and in this section we will analyze the different 
parameters that influence the modulation bandwidth. 

In Fig. [4] b) we show the modulation bandwidth as 
a function of the pumping rate. The numerically exact 
modulation bandwidth obtained using the master equa- 
tion, as described in appendix |Aj is compared to the 
bandwidth found using the transfer function in Eq. ( 15 ) 
which is based on the rate equations. The figure shows 
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the same three regimes as Fig. [2j We see that the rate 
equation model agrees very well with the master equa- 
tion based model, except for a slight underestimation of 
the modulation bandwidth at intermediate values of the 
pumping rate. However, both models show a singular- 
ity in the modulation bandwidth at a certain value of 
the pumping rate in the lasing regime. The position of 
this singularity corresponds to the pumping rate at which 
the cavity population is at its maximum (see fig. [2]). At 
this point a small increase of the pumping rate, will not 
change the population in the cavity since the slope of n a 
is zero. Thus, the output modulation amplitude is zero, 
and the bandwidth is ill-defined. It should be noted that 
this specific pumping rate also indicates the point where 
an increase of the pumping rate goes from giving a posi- 
tive change of the cavity population to a negative change, 
as evidenced by the slope of the cavity population. 

The constant 7p in the expression for the modulation 
transfer function, Eq. (151, describes the variation of 



F when the pumping rate is changed. This variation is 
small and 7p = can be assumed with only a minor 
change to the result. The modulation transfer function 
then becomes 



H(u) 



(17) 



where u> R = laalee+lealae is the determinant and j m = 
lee + 7oo is the trace of the matrix in Eq. ( fl4| ) . Since we 
have ignored the pump dependence on F, this expression 
is similar to that of a conventional laser [44]. Setting 
|-ff(^3dB)| 2 — 3> corresponding to a 3 dB decrease of the 
response, gives the modulation bandwidth 
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FIG. 4: (color online) a) The probability that the excited state of 
the emitter is occupied, Tig . cLS £L function of the pumping rate. The 
solution using the master equation (solid line) and the rate equa- 
tions (dashed line ) is shown, as well as the analytical approximation 
(marks), Eq. {24} , which is valid in the lasing regime, b) The mod- 
ulation bandwidth as a function of the pumping rate. The modu- 
lation bandwidth is found using the modulation response H(lo) as 
deduced from the master equation results using the method out- 
lined in appendix [A] (full line) or from the rate equations model 
using Eq. | |15| l (dashed line). Also shown is the analytical expres- 
sion from Eq. \\9\ (marks) . The parameters are g = 300 ns _1 , 
K = 100 ns" 1 and 7 = 10 ns" 1 . 



/3dB 



1 

~ 2^ 

V 

2 im 



(18) 



In the limit where 7^ 3> ui% a 



where top = ujf 
second order Taylor expansion of Eq. ( 18 ) in the param- 
eter gives 



7; 



/3dB 



Wfl 1_ laalee + 7ea7a 

2-7T 7 m 2-7T 7 ee + -f aa 



(19) 



which is the final analytical expression for the modula- 
tion bandwidth. The discriminant in Eq. ( 16 ) can be 
written D = 7^ — 4w%, and thus the criterion 7^ ^> bj 2 R 
corresponds to an over-damped system. In this limit, the 
poles of the characteristic equation are 



si 



UJ-r 



7" 



and s 2 = 7n 



In 



In 



(20) 



where we see that s\ exactly corresponds to the mod- 
ulation bandwidth found above. This shows that the 
modulation bandwidth of the system is determined by 
the lower of the two poles of the system. These poles 
are different from those of a conventional semiconductor 



laser, where the poles normally form a conjugated pair 
of complex numbers |44j . By analyzing the modulation 
bandwidth it is thus possible to gain insight into the fun- 
damental dynamic properties of the nanolaser. 

In Fig. [4] b) the modulation bandwidth obtained with 
the approximate expression in Eq. ( 19 ) is compared with 



the numerically exact result found using the master equa- 
tion model, as well as modulation bandwidth found us- 
ing Eq. (15). It can be seen that the two rate equation 



based results are very similar, justifying the approxima- 
tion of disregarding jp. Using the approximate formula 
we can analyze the modulation bandwidth in the differ- 
ent regimes; low, intermediate and high pumping rate. 



A. Low pumping rate 

In the limit of low pumping rate (regime I in Fig. EJ), 
we approximate n a sa n e ~ 0, and the modulation band- 
width becomes 



/3dB 



1 4# 2 (k + 7) + n-f(n + 7) 
2^ 8g 2 + (k + 7 ) 2 ' 



(21) 



using Eq. (19 1. In the limit where k — > and g 3> 
7 we have /3dB j^7- Likewise, when 7 — > and 
j » k we have /3dB — > j^k. And for k — > 00 (7 — > 
00) we have / 3dB -> ^ (?3dB 5^) as P -> 0. This 
shows that an excess cavity population Sn a can decay 
through two channels; either directly out of the cavity 
at rate k, or when the emitter and cavity are strongly 
coupled, it can reexcite the emitter which can then lose 
the excitation trough spontaneous emission at the rate 
7. In this respect, the emitter effectively functions as a 
drain for the cavity population. 



B. Lasing 

When the system is lasing (regime II in Fig. [4]) the 
emitter is inverted, meaning that the probability of the 
emitter being in the excited state is larger than the prob- 
ability of it being in the ground state, i.e. n e > n g . In 
Fig. [4] b) we see that the modulation bandwidth de- 
creases as the emitter becomes increasingly inverted. We 
can examine this using Eq. (19 1. "f m does not change 



much in this regime, so we can concentrate on the deter- 
minant uj 2 r . Inserting j aa and j ee and using Eq. (13 1 the 
determinant can be written 



OJr 



F(P-f + K + T(n g -n e ) 



kT. 



(22) 



The last term inside the parentheses becomes negative 
when the emitter becomes inverted, and is thus responsi- 
ble for the decrease in the modulation bandwidth. This 
effect would not be seen in conventional semiconductor 
lasers where, as discussed above, T is independent of P 
and small compared to the other rates. Thus we again 
see that the inclusion of the pumping rate in T has a large 
impact on the dynamics of the system. 



The solution to the rate equations Eq. (10) - (11) can 



be found analytically in the limit of small n. The cavity 
population is described by [T5] 



P-7 (P + 7) 5 



2k 8g 2 
and the excited state occupation is 

1 k(P + 7 ) 



(23) 



(24) 



In the lasing regime the probability of the emitters being 
in the excited state, grows linearly from the 
pumping rate is increased. This approximation is shown 
in Fig. [t] a). Inserting Eq. Q and ([24) into Eq. ([Ej} 
the modulation bandwidth becomes 



/3dB 



K 4g 2 (P - 7 + k) + kT{3k - T) 
2tt 4g 2 (P - 7 + k) + 2k 2 (k + T) '' 



(25) 



where it is apparent that the modulation bandwidth is 
very dependent on the cavity decay rate. However, since 
7 <C P in this regime, the modulation bandwidth is al- 
most independent of the background decay rate. 



C. High pumping rate 

For very large pumping rates, corresponding to regime 
III in Fig. [4j the emitter is fully inverted (n e = 1), but 
because of self-quenching the cavity is empty (n a = 0). 
The modulation bandwidth is therefore 



/3dB 



1 F(k-T) + kT 



2tt 



k + t 



(26) 



which tends to /3dB ~ k 5 - when P is large. We see that 
the expression does not depend on 7. Unlike in the low 
pump rate limit, the emitter is fully inverted and cannot 
absorb a photon from the cavity, thus the emitter cannot 
function as a decay path for the excess cavity population. 
This cavity population can therefore only decay out of 
the cavity at rate k to bring the system back into steady 
state. 



V. MODULATION BANDWIDTH FOR 
SYSTEMS WITH MANY EMITTERS 

So far, we have calculated the modulation response 
only for a system with a single emitter. However, it 
is straight forwards to adapt the rate equation based 
expression for the modulation response to take several 
emitters into account. We see from Eq. ( [l0| that if 
the emitters are similar, the only change in the equation 
for n a when going from one to N emitters is F —¥ NF. 
This corresponds to a change of the light-matter cou- 
pling strength from g to y/N g as described previously. It 
means that in the case of TV-emitters, the elements of the 
matrix for the linearized rate equations are 



7ae = iVP(2n a + l), 

7ea = P(2n e -1), 

7aa = Nj ea - K, 

lee = -lae/N - 7, -P 



(27) 
(28) 
(29) 
(30) 



( 19 ) to obtain the 



which can be used directly in Eq. 
modulation bandwidth. 

In Fig. [5] the modulation bandwidths as a function of 
the pumping rate per emitter for up to 25 emitters are 
shown. These results have been obtained using the rate 
equation based method generalized to TV-emitters. In 
general we see that the form of the curves is similar to the 
one emitter case. The systems with many emitters show 
a larger dip in the modulation bandwidth, and in addi- 
tion the pumping rate at which the emission quenches is 
larger, making the lasing regime broader. Also for more 
emitters, there seems to be a broader plateau of high 
modulation bandwidth before the dip, which indicates 
that these systems are less susceptible to the bandwidth 
lowering effect of the pumping method. However, the 
dips are deeper when the number of emitters is increased. 
The figure also shows that using more emitters does not 
increase the largest possible modulation frequency, since 
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Pumping rate per emitter, P / g 



FIG. 5: (color online) The modulation bandwidth as a function of 
the pumping rate per emitter for N = {1, 2, 5, 10, 25} obtained us- 
ing the rate equation model generalized to TV-emitters. The param- 
eters are gi = g = 300 ns _1 , k = 100 ns — 1 and 7; = 7 = 10 ns — 1 



same behavior as the single emitter case, although the 
system can operate at higher pumping rates and sustain 
higher power levels, before the output quenches or the 
dip in the modulation bandwidth occurs. 

In this work we have modeled the pumping process as 
an incoherent rate from the ground to the excited level of 
the emitter. We have seen that the pumping process has 
a large influence on the dynamics, and that because the 
pumping induces a large dephasing of the photon assisted 
polarization, it is possible to write a set of rate equation 
for the system. The same dephasing was found to be 
responsible for the quenching of the output. We have also 
shown that the pumping rate enters the rate equations 
in a different way than in conventional semiconductor 
lasers, giving an additional decay of the excited state of 
the emitter. This extra term in the rate equations is 
responsible for a dip in the modulation bandwidth and 
the over-damping of the response. 



it is limited by the cavity decay rate k. However, with 
more emitters it is possible to obtain lasing in a system 
with a small cavity Q-value, i.e. larger k, giving a corre- 
sponding larger modulation bandwidth. 

A generalization of the results in section |IV A| and 



IV C 



shows that /3dB 



JV+127T' (/3dB — * Tv+l 27i 

(k — > 0) in the low pumping rate regime. In Fig. [5] 



2i) aS 



7 ' 



we recognise the 



^y3Tj behavior which decreases the mod- 
ulation bandwidth in systems with more emitters. In the 
high pumping rate regime we have /3dB — > which is 
similar to the one emitter case. 



2tt : 



VI. SUMMARY 

In this article, we have derived a set of rate equations, 
which can be used to model the behavior of nanolasers 
based on single or few two-level emitters. The steady 
state values of the cavity population and emitter occu- 
pation, have been compared with result obtained using a 
full master equation model. The models agree when the 
systems are in the lasing regime and when the systems 
are quenched at large pumping rates. 

Using these rate equations, we have calculated the dy- 
namical properties of the nanolasers. We have shown 
that, contrary to conventional semiconductor lasers, 
these nanolasers typically respond in an over-damped 
way when subject to a modulation of the pumping rate. 
Using the rate equations, we have also been able to ob- 
tain analytical expression for the modulation bandwidth 
of the lasers, and shown they exhibit a non-trivial behav- 
ior, with e.g. a dip in the bandwidth for a pumping rate 
in the lasing regime. For comparison, we have calculated 
the modulation response and bandwidth from the more 
fundamental master equation using first order perturba- 
tion theory. These results are in good agreement with 
those obtained from the rate equations. 

For systems with several emitters, we have seen the 
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Appendix A: Calculating modulation response using 
perturbation theory 

This appendix outlines how to calculate the modula- 
tion response directly from the fundamental master equa- 
tion using first order perturbation theory. 

The master equation, Eq. (JlJ, for the density matrix 
is linear and can thus be written as 



dp 
~dt 



(Al) 



where p is a vector containing the elements of the den- 
sity matrix. The steady state solution, p ss , is given by 
Mp ss = 

When the pumping rate is time-dependent the system 
matrix, M, becomes time-dependent. If the system is 
perturbed by a small change in the incoherent pumping 
rate, P, the matrix can be written as M + oMi(t), where 
M is time- independent and aMi(t) specifies the small 
perturbation. Using a series expansion for p 



pit) = 5>y ">(*), 



71=0 



the zero order solution is given by 



dt 



= M3(0) 



and the first order solution can be found from 



dt 



(A2) 



(A3) 



(A4) 
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By writing (f t ~M)^(t) = Mi(i)p< )(i) it is apparent 
that Mi (t)p^ ^ (t) is the driving term. Defining the Greens 
function as 



d_ 
dt 



(A5) 



G(t-i') = I5(t-t'), 
the first order correction can be written 

fl l \t) = /eft'G(t-*OM 1 (t')p (0) (*')- (A6) 



The Greens function is found by Fourier transforming 
both sides of ( A5 ) , and is 



G(w) 



1 



(A7) 



Identifying the right side of Eq. ( A6 1 as a convolution, 
the Fourier transform of p^ 1 ' is 

tfU (w) - G(w)F t (Mx (t)^ (i)) , (A8) 

where F t is the operator for the Fourier transform with 
respect to t. The cavity occupation is given by 

n cav (t) = Tr [a^ap(t)] 

= (9,n\p(t)\g,n) + (e,n\p(t)\e,n) ) 



n,\G{e,g} 



(A9) 



Thus the Fourier transform of the first order correction 
to the cavity occupation is 



>»= n Pn!x^)' (A10) 

n,\G{e,g} 



The matrix Mi(t) 



XiJ(t — i ss ), where t ss is 



the point in time after which the system is assumed 
to be in steady state, specifies the perturbation, and 
contains a delta function in the elements where M 
contains the pumping rate P. All other elements 
are zero. Thus the Fourier transform in (A8) be- 

comes X x F t (^(t)S(t-t. B )) = X^ ) (t„) exp{-iut ss ). 
Therefore, the modulation transfer function, defined as 
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